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6 ABSTRACT 

+-» , The evolutionary properties and spatial distribution of I Zwicky 18 stellar populations are analyzed 

by means of HST/ACS deep and accurate photometry. The comparison of the resulting Colour- 
Magnitude diagrams with stellar evolution models indicates that stars of all ages are present in all 
the system components, including objects possibly up to 13 Gyr old, intermediate age stars and very 
young ones. The Colour-Magnitude diagrams show evidence of thermally pulsing asymptotic giant 
branch and carbon stars. Classical and ultra-long period Cepheids, as well as long period variables 
have been measured. About 20 objects could be unresolved star clusters, and are mostly concentrated 
in the North- West (NW) portion of the Main Body (MB). If interpreted with simple stellar population 
' models, these objects indicate a particularly active star formation over the past hundred Myr in IZw 18. 

The stellar spatial distribution shows that the younger ones are more centrally concentrated, while 
\Q . old and intermediate age stars are distributed homogeneously over the two bodies, although more 

easily detectable at the system periphery. The oldest stars are best visible in the Secondary Body 
(SB) and in the South East (SE) portion of the MB, where crowding is less severe, but are present 
also in the rest of the MB, although measured with larger uncertainties. The youngest stars are a few 
Myr old, are located predominantly in the MB and mostly concentrated in its NW portion. The SE 
portion of the MB appears to be in a similar, but not as young evolutionary stage as the NW, while 
the SB stars are older than at least 10 Myr. There is then a sequence of decreasing age of the younger 
r ^ ' stars from the Secondary Body to the SE portion of the MB to the NW portion. 

' All our results suggest that IZwl8 is not atypical compared to other BCDs. 

Subject headings: galaxies: evolution — galaxies: stellar content — galaxies: dwarfs — galaxies: 
individual (IZwl8) 



1. INTRODUCTION 

I Zwicky 18 (also Mrk 116 or UGCA 166, hereafter 
IZwl8) is the most intriguing Blue Compact Dwarf 
(BCD) galaxy in the local Uni verse. It became f amous 
right after its discovery, when ISearle fc Sargenti (|1972fi 
measured from its emission-line spectrum an oxygen 
abundance [0/H]= -1.14 dex, only 7% of the solar value, 
indicating an almost unprocessed gas content. The first 
studies on its color and composition (|Sargent fe Searld 
[19701: ISearle fc Sargenti H97l ISearle et al.lll973D empha- 
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sized its almost primordial evolutionary status. IZwl8 
shows very blue colo rs, U— B= —0.88 and B— V= —0.03 
(jvan Zee et al.l fl998h . suggesting the presence of a very 
young stellar population, with a current star forma- 
tion rate (SFR) much higher than the past mean value. 
The total mass of IZwl8 from the rotation curve at 
a rad ius of 10"-^12" is estima t ed to be ~ 10 * * M^, 
(e.g.. iDavidson fc Kinmanl 119851 : iPetrosian et al.1 119971 : 
Ivan Zee et al.lll998f ). where the neutral gas corresponds 
to ~70% of the total mass, but only 10 7 M Q of HI 
is associated with the optical part of the gala xy (e.g., 
iLequeux fc Viallefondl[T98Cl : Ivan Zee et al.lll998fl . 
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metal deficiency around 1/30 
IZwl8 remains the star forming galaxy with the sec 
ond lo west metallicity an d the lowest helium content 
known (jlzotov et al.l 120091 ) . This makes the system a 
fundamental point in the derivation of the primordial 
helium abundance (llzotov et al.lll994t IQlive et al.lll997l; 
llzotov fc Thuanlll998l ) and in the study of the properties 
of chemically unevolved galaxies. 

When discovered by Zwicky (1966), IZwl8 was de- 
scribed as "two galaxies separated by 5". 6 and inter- 
connected by a narrow luminous bridge" , surrounded 
by two "very faint flares" at 24 " northwest. Subse- 
quent CCD ground-ba sed images (jDavidson et a,l.|[l989t 
iDufour fc Hester! I199O0 revealed a more complex struc- 
ture: the two galaxies turned out to be two star- 
forming regions of the same galaxy (usually indicated 
as NW and SE components), while the two flares are 
just the most prominent of a few nebulosities surround- 
ing IZwl8. These minor systems are roughly aligned 
toward the northwest and were initially considered at 
the same distance, but subsequent spectroscopic stud- 
ies showed that only one component (referred to as 
component C in IDavidson et al.l I1989D is at the same 
distance as IZwl 8 and is physically associated with 
the main body (IDufour et al.l Il996at iPetrosian et al.l 
119971; Ivan Zee et al.l 11998ft . The other diffuse objects 
have been recognize d as background galaxies (see, e.g., 
IDufour et aT1ll996bH . In the following we will refer to 
IZwl8 and to component C respectively as main body 
(MB) and secondary body (SB) of IZwl8. Both bodies 
have been resolved into single stars for the first time by 
H ST with the Wide Fie l d Pla neta ry Camera 2 (WFPC2 ) 
bv lHunter fc Thronsonl (|1995l ) and fDufour et al.l (jl996b[ ). 

The observational features described above triggered 
the key question over the nature of IZwl8: is it a young 
galaxy which is presently experiencing its first burst of 
star formation, or is it an old system which has already 
formed stars in the past? HST studies by different au- 
thors have tried to characterize the evolutionary status of 
IZwl8, with rather discordant results, but with a clear 
trend intriguingly reminiscent of Oscar Wilde's Dorian 
Gray's pi cture: the more o ne looks at it, the older it 
appears. IMoisi et al.l (|1999h . using HST/WFPC2 data, 
were the first to detect Asymptotic Giant Branch (AGB) 
stars in the galaxy, and thus to demonstrate that IZwl8 
is as old as at least several hundre d Myr, and possibly 
much older, llzotov fc Thuanl (|2004f ) failed to detect red 
giant branch (RGB) stars with HST observations made 
by the Advanced Camera for Surveys (ACS), and thus 
concluded that the galaxy is at m ost 500 Myr old. How - 
ever, both iMomanv et al.l (|2005D and lTosi et all (|2007l ). 
from independent re-analyses of the same ACS data set, 
suggested that IZwl8 should be older than at least 1-2 
Gyr, since it did appear to contain also Red Giant Branch 
(RGB) stars. O ur own time-series HST/ACS photome- 
try (GO 10586. [MoTsTet~aTl [20071 hereafter A107), has 
allowed us to shed light on the situation, thanks to the 
combined information from the color-magnitude diagram 



(CMD) and the lightcurves of the Cepheids we were able 
to measure. 

In fact, A107 not only characterized the RGB and iden- 
tified its Tip (TRGB) at I = 27.27±0.14mag, but, more 
importantly, detected for the first time a few Cepheids 
whose lightcurves allowed to independently and more re- 
liably pin down the galaxy's distance to the couple Mpc 
accuracy (Fiorentino et al. 2010, hereafter F10; Mar- 
coni et al. 2010). The distance to IZwl8 was estimated 
from the Cepheids using theoretical V, I Wesenheit rela- 
tions specifically computed for the very low metal abun- 
dance of the galaxy, and by theoretically fitting the ob- 
served lightcurve of the most reliable Classical Cepheid 
with period 8.71 days. The Wesenheit relations provide 
an intrinsic distance modulus for IZwl8 of (m — M)q = 
31.4±0.2 (D=19.0l^'7 Mpc) adopting a canonical mass - 

luminosity relation, and of 31.2 ±0.2 (0=17.41^ Mpc) 
when an over luminosity of 0.25 dex is assumed for each 
given mass (non canonical scenario). The theoretical 
modeling of the 8.71 days Cepheid lightcurve provides 
31.4 ±0.1 (D=19.0±^ Mpc). The TRGB identified by 
A107 implies an instrinsic distance modulus of 31.30 ± 
0.17, in excellent agreement with the Cepheid-based es- 
timates. This agreement, in turn, provides independent 
support to our interpretation of the CMD faint red se- 
quence as an RGB. 

Hence, IZwl8 is a ctually farther aw ay than originally 
thought (for instance lAloisi et aT]|1999l adopted from the 
literature a mean value of D=10 Mpc) and is in prac- 
tice at slightly larger distance than the Virgo cluster of 
galaxies. This made it impossible in the past to resolve 
its older, hence fainter, stars, and quite challenging even 
with the outstanding performances of HST's current in- 
strumentation. Today IZwl8 is one of the most distant 
system ever resolved into stars. 

In this work, we have used the photometric catalogues 
of our previous works, which have been focused on de- 
termining the age (A107) and the distance through the 
variable stars (F10). In this paper we describe the evo- 
lutionary properties and the spatial distribution of the 
resolved stellar populations and candidate star clusters 
of the main and secondary bodies of IZwl8. In Section 2, 
we briefly summarize the main differences in the photo- 
metric procedures adopted by A107 and F10. The CMDs 
and the spatial distribution of stars in different evolution- 
ary phases are presented in Section 3 and 4, respectively. 
The results are analysed in Section 5 and summarized in 
Section 6. 

2. THE PHOTOMETRIC DATA-SET 

The data acquisition and reduction have been exten- 
sively described by A107 and F10. Here we only sum- 
marize their main aspects. Time series photometry of 
IZwl8 in the F606W (broad V) and F814W (I) filters 
was obtained with the Wide Field Channel (WFC) of 
the ACS in 13 different epochs properly spread over a 
time interval of 96 days (program Id. 10586, PI Aloisi). 
These proprietary data were complemented by 5 F555W 
(V) and 3 F814W archival ACS/WFC images of IZwl8 
from HST (program Id. 9400, PI Thuan). 

IZwl8 is quite distant, plenty of bright stars hiding 
and/or contaminating the fainter ones, and plenty of gas, 
which complicates the photometry of individual stars. 
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In addition, the throughput of the F606W filter is sig- 
nificantly affected by the H a contribution of ionized gas 
(quite abundant in the MB of IZwl8) and makes the star 
detection and measurement more complicated. Thus, de- 
pending on the goal, one can find more appropriate to 
restrict the analysis of the galaxy populations to safer 
but fewer objects (A107), or to more uncertain but nu- 
merous ones (F10). 

In A107, to produce the galaxy CMD, we built four 
coadded images by stacking all single epoch images per 
filter for both archival and proprietary data, separately. 
Then we performed PSF photometry of each of these 
four coa dded images u sing the DAOPHOT / ALLSTAR 
package (|Stetsonlll987f) . The calibration to the Johnson- 
Cousins photometric syst em was obta i ned u sing the 
transformations given by iSirianni et al.l (|2005[ ). Since 
the goal of that paper was to demonstrate that faint red 
stars actually existed where none had previously been 
detected, the published CMD was produced following a 
strategy aimed at minimizing the number of false detec- 
tions. As a result, it contained only the 2099 sources 
which were independently detected in all the four deep 
images (the Vs and Is of both datasets) , with good match 
in the position (within one pixel) and magnitude (within 
3.5cr) of each source in the different images, and with 
sharpness \sharp\ <0.5. These are the most reliable ob- 
jects measured in each of IZwl8's bodies. 

In the case of F10, to produce the lightcurves 

of the candidate variables, we made u se of 

DAOPHOTII/ ALLSTAR/ ALLFRAME (jStetsonl H987L 
1994) to get PSF photometry on each single epoch (15 
/ and 18 V images, in total between proprietary and 
archival data). To study the variability we relaxed the 
A107 selection criteria, since variable stars are more 
easily distinguishable from spurious objects thanks to 
their brightness variation in the images at different 
epochs. Moreover, to collect a good sample of candidate 
variables at all magnitudes, we examined objects close 
to the detection threshold. This photometry resulted in 
a J, V — I CMD of IZwl8 which contains 6859 sources 
with photometric errors av,i < 0.5 mag and sharp 
parameter in the range from — 1 to + 1 , and reaches V ~ 
30.0, i.e. more than half a magnitude deeper than A107's 
photometry. The calibration to the Johnson-Cousins 
photometric sy stem was obtai n ed us ing the transforma- 
tions given by ISirianni et al.l (|2005T) . as for the other 
data set. 

As shown above, there is a remarkable difference in the 
number of objects selected with the two approaches. To 
infer the properties of the stellar populations inhabiting 
IZwl8 minimizing the effects of uncertain photometry, 
we will base our results on the hyper-selected A107 cat- 
alogue. Nonetheless, we consider it useful to show some 
features of the larger F10 catalogue both for complete- 
ness and for comparison. 

Since our aim is to analyse the evolutionary properties 
of IZwl8's stellar populations, we have put special care 
in verifying whether the selection differences affect some 
CMD regions more than others. We have hence divided 
the CMDs of the two sets in 8 zones roughly represen- 
tative of different evolutionary phases. Since the TRGB 
was detected by A107 at / = 27.3, we have adopted this 
magnitude as the border between bright and faint stars: 
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Figure 1. CMD of the A107 hyper-selected stars of the main 
(upper panel) and secondary (lower panel) body showing the 8 
evolutionary zones described in the text. Notice that zones 3 and 
5 partially overlap: this is intentional to let the brightest stars of 
the blue plume be considered either as upper MS objects or as 
supergiants or as star clusters. 



• zone 1, stars in the extremely hot portion of the 
blue plume with V — I < —0.2 and brighter than 
I = 27.3; 

• zone 2: stars in the extremely hot portion of the 
blue plume with V — I < —0.2 and fainter than 
/ = 27.3; 

• zone 3: stars in the "standard" blue plume, with 
-0.2 < V - I < +0.2, and brighter than I = 27.3; 

• zone 4: stars in the "standard" blue plume, with 
-0.2 < V - I < +0.2, and fainter than I = 27.3; 

• zone 5: objects with I < 23 and any colour; 

• zone 6: stars with V-I > 0.8 and 24.5 < I < 27.3; 

• zone 7: stars with V — I > 0.8 and / > 27.3; 

• zone 8: all the remaining ones;. 
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Figure 2. Portion of one F814W image of IZwl8, with overimposcd the 6859 stars measured by F10. The X-Y coordinates are in pixels 
(pixel scale ~ 0.035" /pix). The white cross indicates the center of this portion ACS/WFC image, which has RA=09 h 34 m 01 s , DEC=55° 
14' 33.3". The stars are divided and color-coded according the 8 CMD evolutionary zones in Fig.[T] 



The number of objects measured in each zone is reported 
in Table 1 for both bodies and both data sets. 

The zones are displayed on the A107 CMDs in Figure[TJ 
Figure [2] shows an F814W image of IZwl8 with the F10 
measured objects colour-coded according to the 8 CMD 
zones of Figure [1] We arbitrarily define as MB objects 
those with Y coordinate in the image smaller than 1200 
px and as SB objects those with Y>1200 px. In FigureQ] 
and all the following ones, the CMDs of the two bodies 
are splitted following this definition. 

We have found that the two photometries agree well 
with each other, once we take the different selection cri- 
teria and photometric depth into account. The objects 
measured by both A107 and F10 fall consistently in the 
same CMD zone (i.e. have the same colours and magni- 
tudes). 

The F10 objects not present in the A107 hyper-selected 



catalogue fall mostly in the fainter portions of the CMDs, 
as expected since they can have larger photometric error. 
An intriguing difference is that the MB zones 1 and 2 
are much more populated in F10 than in A107 (see Ta- 
ble 1), with a ratio N(F10)/N(A107)=8.5 in zone 1 and 
10.9 in zone 2, significantly higher than the total MB 
ratio 4744/1187=4.0. Figure [5] reveals that most of the 
F10 objects of zones 1, 2 and 4 fall on the gas-rich regions 
in and around the MB. A possible cause for the discrep- 
ancy is therefore gas contamination in the F606W flux, 
even if in principle all photometric packages should over- 
come this kind of effects by subtracting the surrounding 
background flux to each source. 

As a matter of fact, many of the objects in zone 1 
and 2 show a discrepancy in ^-magnitude depending on 
the used filter, e.g. F606W or F555W. If, we restrict 
the selection criteria on the F10 catalogue and keep only 
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Table 1 

Number of MB and SB objects in the A107 and FIO data sets, and 
corresponding percentages, in the 8 CMD zones defined in the text 
and shown in Figure [T] 



CMD 
zone 


Na/s 
(A107) 


% 


Na/s 
(FIO) 


% 


(A107) 


% 


N S s 
(F10) 


% 


zone 1 


76 


6 


649 


13 














zone 2 


84 


7 


917 


20 


43 


5 


108 


5 


zone 3 


277 


23 


820 


17 


93 


9 


98 


4 


zone 4 


186 


16 


707 


15 


262 


29 


584 


28 


zone 5 


34 


3 


51 


1 


2 


<1 


2 


<1 


zone 6 


105 


8 


210 


1 


60 


6 


110 


5 


zone 7 


126 


11 


349 


8 


120 


13 


375 


18 


zone 8 


305 


26 


1060 


22 


332 


36 


838 


40 


total 


1187 


100 


4744 


100 


912 


100 


2115 


100 



a Noticc that the total number of stars in the MB is smaller than the 
sum of the values in the corresponding column, because the objects 
falling in the overlapping zones 3 and 5 are counted twice:19 in F10 
and 6 in A107. 

the objects with ay.i <0.1, most of those in the CMD 
zone 2 and several of those in zone 1 are rejected. We 
thus consider as real stars on the gas filaments only the 
brighter ones with smaller photometric error. 

Finally, we have individually and visually inspected the 
brighter objects of FIO zone 1 using appropriate routines, 
and found that some of them do appear as real stars. We 
consider the vast majority of the blue objects measured 
by F10 on the gas filaments and not present in the A107 
catalogue as likely spurious detection due to gas peaks, 
but we do believe that some of the brighter zone 1 de- 
tections are real stellar objects. 

To be conservative, we will base all our results on the 
restricted A107 catalogue. 

3. COLOUR-MAGNITUDE DIAGRAMS 

To let the reader easily figure out the involved stel- 
lar ages, the left panel of Figure El shows the CMDs of 
the hyper-selected A107 objects with o verimposed theo- 
retical isochrones (|Girardi et al.l [2010h . The isochrones 
have metallicity Z=0.0004 and are plotted adopting 
TO - M)n=31.4 (A 107 and FIO) and E(B-V)=0.032 
van Zee et all 1 19981) . The CMDs corresponding to the 
less restrictively selected F10 photometry are displayed 
in the right panel of Figure El where stellar evolu- 
tion tracks are also plotted to let the reader visualize 
what the involved star masses are. The plotted tracks, 
ag ain with Z=0. 004, are the classical Padova models 
bv lFaeotto et al.l (|1994D . We note that IZwl8 is located 
at high galactic latitude (b = +45°) and therefore its 
CMDs do not suffer from significant contamination from 
foreground stars belonging to our Galaxy, meaning that 
is not necessary to correct for it our observed diagrams. 

IZwl8 is a starforming dwarf and its CMDs show all 
the features typical of these galaxies: a prominent blue 
plume with mean V — I ~ — 0.1 and a fainter and broad 
red plume with V — I > 0.8. Once compared with the 
stellar evolution tracks of Figure El the blue plume turns 
out to correspond to high- and intermediate-mass stars in 
the main sequence (MS) phase or at the blue edge of the 
blue loops covered during central He burning. These two 
phases are resolved by HST photometry in nearby late- 
type galaxies, but at the distance of IZwl8 the photo- 



metric error inevitably is large enough to let them merge 
in just one broad sequence. Stars with colours inter- 
mediate between the blue and the red plume are either 
supergiants (brighter than 7 ~ 23) or intermediate-mass 
stars in the blue loop phase. The brighter ones might 
also be unresolved star clusters. 

The red plume is populated by red supergiants (7 < 
23), AGB (24.5 < 7 < 27.3) and RGB stars (7 > 27.3). 
Our data allow to cover the brighter 1.5 magnitude por- 
tion of the RGB, where UDaophot is ±0.1 mag in 7 and 
±0.2 in V — I (A107) and the photome tric error, esti- 
mate d from extensive artificial star tests (jAnnibali et al.l 
1201 If) , is smaller than ±0.3 mag in 7 and ±0.4 in V — I 
(see Figure []}. The colour extension of the RGB is larger 
than this error and reflects, as usual, both the spreads 
in age and in metallicity of the resolved stars. However, 
given the extremely low metallicity of IZwl8, in this case 
it is mostly the age that contributes to the RGB thick- 
ness. In other words, the reddest RGB stars are most 
likely low mass stars 10-13 Gyr old. 

The redder red plume stars with 7 ~ 26.4 are presum- 
ably carbon stars, not too differe nt from those measured 
in th e Magellanic Clouds by Ivan der M arcl & Cioni 
(|2001| ). The redder objects at 7 ~ 25.5 are most proba- 
bly stars of 4-5 M Q in their thermal pulsing (TP) phase 
at the end of the AGB, since their location in the CMD 
is perfectly consistent with the predictions of Marigo & 
Girardi's (2007) TP-AGB models. These authors show 
in fact that at very low metallicities (Z < 0.001) in (rela- 
tively) massive intermediate-mass stars efficient nuclear 
burning at the base of the convective envelope may favour 
the formation of C-rich models and make their tracks 
describe pronounced excursions towards lower effective 
temperatures. These excursions are recognizable in Fig- 
ure [3] in the tangled, messy portions of the isochrones 
above the TRGB. 

Not only does the MB, by definition, host a larger num- 
ber of stars; it also contains a larger fraction of young, 
blue and bright stars than the SB. Both the luminosity 
and the number of the brighter stars in the SB are sig- 
nificantly lower than in the MB. Almost no supergiant 
is found in the SB and no Cepheid. Moreover, while the 
blue plume of the MB covers the entire range of blue mas- 
sive stellar evolution tracks and youngest isochrones, that 
of the SB turns off the MS around 7 ~ 28, with appar- 
ently no bright star younger than 10 Myr being detected 
there. This is a sig nificant difference, which reflects, as 
discussed already bv lAloisi et al.l (|1999| ). both a different 
rate and a different history of star formation, with the 
MB being more active than the other at recent epochs, 
and the SB more active earlier than about 100 Myr ago, 
as suggested by the higher density of objects around the 
isochrones older than this age in the left-bottom panel 
of Figure El At the faintest magnitudes, we should also 
consider that crowding affects the MB more severely, pre- 
venting the detection of its faintest stars. This may be 
particularly the case in the NW part of the MB. 

The variable stars in IZwl8 have been identified by 
F10 exploiting our ACS time-series photometry. The 
intensity-averaged mean magnitudes and colors of the 
confirmed Cepheids are plotted as filled (red) triangles in 
Figure El while those of the confirmed Long Period Vari- 
ables (LPVs) are plotted as filled (red) squares. Their 
position (see F10 for identifications and nomenclature) 
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Figure 3. CMD of IZwl8 MB (top panels) and SB (bottom panels) obtained by A107 (left panels) with restrictive selection criteria and 
by F10 (right panels) with less restrictive ones (see text). The CMDs are calibrated in the Johnson-Cousins photometric system. The 
three confirmed Cepheids are highlighted by filled (red) triangles, and the two Long Period Variables by filled (red) squares. All confirmed 
variables are plotted according to their intensity-averaged mean magnitudes and colors from F10. In all panels, stellar evolutionary models 
with metallicity Z=0.0004 are overi mposed to the dat a adopting (m — M)o=31.4 and E(B-V)=0.032. In the left panels we plot isochroncs 
with detailed TP-AGB modelling HGirardi c t al. 2010), and log(age) 6.75, 7.00, 7.25, 7.50 ... 9.75, 10.00 (from top to bottom). The tangled, 
messy portions of the isochrones above the RGB correspond to the TP-AGB unstable phases. In the right panels we plot the classical 
evolutionary tracks by Fagotto ct al. (1994). In green (in the electronic version) are the tracks for 3, 4, 5, 7, 9, 12, 15, 20, 30 and 60 Mq, 
and in red the 6 Mq one, which best fits the 8.71d Cepheid. The TP-AGB phase is not covered by these tracks. The RGB portion of both 
the isochrones and evolutionary tracks are shown in red in the electronic version. 



in the CMD indicates that the two bright ultra-long pe- 
riod Cepheids (ULPs) VI and V15 are supergiants, un- 
less they are hosted by unresolved star clusters making 
them appear brighter. The Classical Cepheid V6 with 
period P=8.71d is on the blue loo p and is best fitted b y 
a -6 M star (see also F10 and iMarconi et al.l l2010h . 
The brighter LPV V4 is a supergiant, while the other 
confirmed LPV V7 lies at the bright edge of the CMD 
region of the TP-AGB stars and is presumably a Mira. 

For a finer analysis of the evolutionary properties of 
IZwl8 components, we have further divided the A107 
objects according to their position in the MB. Figured] 
shows the CMDs of the objects with Y coordinate in the 
image larger (top panel) or lower (bottom panel) than 
650 px (see Figure [2]). This in practice separates ob- 
jects in the NW (top) part of the MB from those in the 
SE (bottom). The two CMDs are fairly similar to each 
other, although the NW one contains more objects in 
zone 5 and is more dispersed than the other. In particu- 
lar, the red plume of the SE region is much better defined 



than that of the NW and allows to clearly distinguish the 
RGB and its Tip from the AGB. The larger spread of the 
NW CMD is likely due to the presence of a larger num- 
ber of bright objects and thick ionized gas filaments that 
make the measurements of fainter (and older) stars more 
uncertain. It is plausible that the gas ionization is itself 
due to the massive luminous objects populating the NW 
region; so, they may represent the only real difference 
between the northern and southern parts of the MB. We 
will further discuss these objects in the next section. 

In spite of the very different number of objects selected 
in the two data sets, once the different depths of the 
resulting CMDs (Figure [3|) are taken into account, their 
stellar evolutionary phases turn out to be populated in 
similar proportions, both in the MB and in the SB (see 
Table 1). As discussed in the previous section, the only 
significant difference is that the MB blue plume is much 
more populated, bluer and thicker in the F10 CMD. 
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Figure 4. CMD of the MB objects in the hyper-selected A107 cat- 
alogue, separating those in the NW region (Y>650 px, top panel) 
from those in the SE region (Y<650 px, bottom panel). There are 
478 objects in the NW region and 704 in the SE one. 

4. SPATIAL DISTRIBUTION AND EVOLUTIONARY 
PROPERTIES OF IZW18 RESOLVED STARS 

We have analysed the spatial location of the various 
types of stars populating IZwl8 and checked if stars of 
different ages are evenly or unevenly distributed. To 
this purpose, we have used the CMD zones of Figure [TJ 
Thanks to the comparison with the evolutionary tracks, 
we can assign to each zone the following evolutionary 
properties: 

• zone 1 corresponds to massive and intermediate 
mass stars on the MS , younger than ~ 10-15 Myr; 

• zone 2 to fainter intermediate mass stars on the 
MS, younger than~ 70 Myr; 

• zone 3 to massive and intermediate mass stars on 
the MS or at the blue edge of the blue loop phase, 
younger than ~ 70 Myr; 

• zone 4 to fainter intermediate mass stars on the MS 



or at the blue edge of the blue loop phase, younger 
than ~ 300 Myr; 

• zone 5 to supergiants, younger than ~ 10 Myr and, 
possibly, to unresolved star clusters; 

• zone 6 to AGB stars, ages from ~ 35 to 300 Myr; 

• zone 7 to RGB stars, ages from ~ 1 to 13 Gyr; 

• zone 8 in most cases corresponds to intermediate- 
mass stars in the blue loops, but also to some su- 
pergiants, ages from ~ 15 to 300 Myr. 

The spatial distribution in the X-Y map of IZwl8 of 
the A107 objects measured in the two bodies and falling 
in each of the 8 CMD zones is reported in Figure [5] 

This map shows that the SB does not contain stars 
in zone 1. In other words the SB has had very little, if 
any, star formation activity in the last few Myr. Zone 
1 stars covering the whole range of magnitudes (hence 
with masses from 15-20 M Q up) are numerous in the 
MB and show a sort of 8-shaped spatial distribution. If 
we look at the fainter portion of the extremely blue stars 
(zone 2), we find that both bodies contain several of these 
objects and that in the MB they seem to delineate the 
same 8-shape as their brighter counterpart. 

The bright stars in zone 3 have masses above 6 M Q and 
are also present in both bodies. Their spatial distribution 
in Figure [5] forms an irregular ellipsoidal clump of objects 
with significant central concentration (the MB) with 
a comma-shaped companion (the secondary: Zwicky's 
"flare"). This is the overall familiar appearance of IZwl8, 
since these are the stars more easily visible, although not 
resolvable from ground. Once again, the SB stars of zone 
3 don't reach the luminosity of the brightest ones hosted 
in the MB (see Figure [lj : while the whole mass and 
luminosity ranges are covered in the latter, only stars 
less luminous than / ~ 24 exist in the former. Taking 
into account that any Initial Mass Function disfavours 
the most massive stars, we consider it more likely that 
the brighter blue stars in the SB are 12-15 M Q stars at 
the blue edge of the blue loop rather than more massive 
stars. In this case, the star formation activity in the sec- 
ondary must have been quite low, or absent, in the last 
10 Myr. 

Fainter (lower mass) blue plume stars (zone 4) can be 
either MS stars with masses between 6-20 M or stars 
at the blue edge of the blue loop with mass in the range 
3-5 Mq (see the tracks in right panel of Figure [3]). They 
are present in both bodies, and are more abundant in 
the SB than in the MB, where they are not found in 
the central regions and have a spatial distribution com- 
plementary to that of their brighter counterpart of zone 
3. This is clearly due to the severe crowding affecting 
the MB and, particularly, its inner regions, where the 
many bright stars completely hide fainter ones. In the 
SB the recent star formation activity has been less in- 
tense, crowding is less severe, and faint stars are much 
more easily detectable. 

The lower activity of the SB at recent epochs is also 
confirmed by the objects of zone 5. As apparent from 
Figure O we find only 2 such bright objects in the SB 
(one of which is quite far away and may actually be a 
background system), and 31 in the MB. The latter are 
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Figure 5. X-Y map of the A107 hyper-selected objects of Figure [3] splittcd according to the 8 CMD evolutionary zones (see text and 
Fig.[TJ. Objects with Y>1200 are associated to the SB, the others to the MB. In the latter, the NW portion corresponds to Y>650 and 
the SE portion to Y<650. 



fairly confined in the inner regions, and most often in the 
NW portion of the MB. They significantly contribute to 
the severe crowding affecting the central area. If inter- 
preted in terms of individual stars, these objects are su- 
pergiants and confirm that the central region of the MB 
has experienced a strong star forming activity in the last 
10 Myr or so. 

Admittedly, some of these very bright objects might 
be not individual stars but unresolved star clusters. We 
discuss this possibility in Section 14.11 

The objects falling in the CMD zone 6 are (see the 
tracks in Figure [3]) AGB stars with mass lower than 9 



Mq, and age (see the isochrones in the same Figure) in 
principle between 35 Myr and 13 Gyr, but in practice 
peaked at, at least, several hundreds Myr, once the stel- 
lar lifetimes and phase durations are taken into account. 
They include the TP- AGB and carbon stars with 4-5 M Q 
(horizontal finger a t I ~ 25.5 in the CMDs) predicted by 
Mari gcTfe Girardil (|2007f) for these low metallicities and 
the carbon stars at / ~26.4 discussed by A107. Zone 6 
objects are equally present and evenly distributed in both 
bodies, although missed in the most crowded portions of 
the MB. 
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The RGB stars are in zone 7: objects with mass be- 
low 1.8-2.0 M Q and age between 1 and ~ 13 Gyr. The 
RGB is notoriously affected by the age-metallicity de- 
generacy, i.e. the ambiguity that stars with the same 
brightness and colour can be older and metal poorer or 
younger and metal richer. Hence, without independent 
spectroscopic metallicity measurements one cannot as- 
sign the age to RGB stars simply from their position in 
the CMD. For this reason the photometric identification 
of the RGB normally puts only a lower limit of 1-2 Gyr 
(i.e. the minimum time required by low-mass stars to 
reach that phase) to the age of the examined stellar pop- 
ulation. However, in the case of IZwl8, its extremely 
low metallicity somewhat relaxes the effect of the age- 
metallicity degeneracy, since the redder RGBs cannot be 
metal rich objects and therefore must be old. In these 
conditions, the colour extension of IZwl8 RGB, which is 
broader than the photometric error and covers the whole 
range of isochrone ages, can be taken as a signature of the 
presence of stars possibly as old as 13 Gyr. RGB stars 
are faint and therefore their detection is strongly affected 
by crowding and photometric performances. This is why 
only with the ACS has it been finally possible to resolve 
and measure them: not even HST was able to reach the 
RGB, with the less sensitive WFPC2. The red giants 
are evenly distributed in the SB and visible only in the 
outer regions of the MB, due to the excessive crowding 
of its inner parts. It is interesting to notice that the MB 
RGBs are located along a semi-arc in the system periph- 
ery. We ascribe this asymmetry to the presence of thick 
gas filaments in the regions where RGBs are not detected: 
overlapping filaments can clearly hide such faint objects. 
We thus believe that RGB stars can actually be evenly 
distributed also in the MB. 

Finally, zone 8 collects all the remaining objects, with 
a large range of brightness and colour, and different ages 
and evolutionary phases. Those with / >24.5 are stars 
in the 3-10 M Q mass range in the central helium burn- 
ing phase and with ages between 30 and 300 Myr. The 
brighter ones are 10-15 M Q mass stars in any post-MS 
evolutionary phase, with age of the order of 15-30 Myr. 
These mixed objects are evenly distributed in both bod- 
ies, with the only exception of the two most crowded 
NW and SE spots of the MB, where they are missing. 
Interestingly, like the objects of zones 2, 4 and 7, they 
are found in similar numbers in the two bodies. Even 
considering the effect of the different crowding, this sug- 
gests that as far as low and intermediate mass stars are 
concerned, the SB is not significantly less populated than 
the MB. 

As for the variable stars, the two ULPs VI and V15 
are located in the NW portion of the MB, while the clas- 
sical Cepheid V6 and the two LPVs V4 and V7 are in 
the SE portion (see F10 figure 5). Interestingly, all the 
confirmed variables belong to the MB: of the 34 good can- 
didate variables examined by F10, almost half (15) were 
actually detected in, or around, the SB, but for none of 
them satisfactory lightcurves and periods were derived. 
Ground-based observations of the brighter variables are 
being performed to better characterize their pulsation 
and evolutionary properties. In particular, further data 
on the ULPs will allow us to verify if they are the exten- 
sion to higher masses of Classical Cepheids and, in case, 
to calibrate them as extremely useful primary indicators 



of cosmological interest. 

Finally, we emphasize that no star is found in the re- 
gion between the MB and the SB, either in the A107 or 
the F10 catalogues. That region looks completely trans- 
parent, with only background galaxies being detected 
there. This shows that the two bodies are completely 
separated. They are members of a two-body system, but 
have no star connecting each other. If tidal interactions 
are at work, they don't appreciably affect the stellar pop- 
ulations of the two components. 

4.1. Cluster candidates 

At 19 Mpc, ACS/WFC photometry cannot recognize 
round regularly shaped clusters whose sizes fall within 
the Point Spread Function. In this section we consider 
the possibility that zone 5 objects are actually not indi- 
vidual stars but unresolved star clusters, and discuss the 
implications on IZwl8's evolutionary status. 

With an intrinsic distance modulus of (m — M)o=31.4, 
the objects of zone 5 have an absolute bolometric magni- 
tude brighter than -8.4, certainly consistent with typical 
cluster luminosities. Actually, the brightest objects, with 

1 ~ 22, hence Mf, j < —9.4, would b e "populous" clus- 
ters, a ccor ding to the definit ions by lLarsen fc Richtlerl 
lj2000f) and lBillett et all (12002ft . These bigger, more mas- 
sive clusters ar e fairly common in starburst dwarfs, suc h 
as NGC 17 05 (iBillett et alj|2002t lAnnibali et all [20091) , 
NGC 1569 (iHunter et all 120001: lOriglia et all 120011) and 
NGC 4449 (jAnnibali et all 120111 ). and tend to be cen- 
trally concentrated, just like these brightest objects of 
IZwl8. Hence, also if taken as star clusters, the objects 
of zone 5 would confirm a higher star formation activity 
in the MB than in the SB. 

To infer the evolutionary properties of unresolved clus- 
ters, two photometric bands are insufficient. Fortunately 
IZ wl8 was obser v ed also in the B band with the WFPC2 
bv lDufour etaD (|1996bD . We cross-correlated the 31 + 

2 objects in Zone 5 with the WFPC2 (F435 W, F555W) 
photometry presented by[Aloisi et~aT](|l999l) . The cross- 
correlation provides 23 + 2 objects photometrized in 
F439W, F555W and F814W. A visual inspection in our 
images of the objects missed by the cross-correlation 
shows that they are all located in the most crowded 
MB region, and thus their position measurement is 
highly uncertain, also because of blending effects. The 
WFPC2 magnitudes were transformed into Johnson- 
Cousins B, V and I ma g nitud es using the transforma- 
tions by[Holtzman et al] (|1995l ). 

In Fig. [5] we show the B — V versus V — I color-color 
diagram and the V versus V — I color-magnitude dia- 
gram for zone 5 objects. The data are compared with the 
Padova simple stellar populations (SSPs) with metallic- 
ity Z= 0.0004, ages from ~1 Myr to ~13 Gyr , and masses 
of 10 3 , 10 4 , and 1O 5 M IjGirardi et alJl2010j) . Zone 5 ob- 
jects have V — I colors in the range —0.2-1.8, and B — V 
colors in the range —0.2-1.4, and are compatible with 
SSP ages from a few Myr to a Hubble time. The two 

objects with V - 1 ~ 0.4, B — V -0.8, and V - 1 0.2, 

B — V —0.4, are located in the most crowded region of 
the main body, and therefore their colors suffer larger un- 
certanties. The objects concentrated around V — I — 0.2 
can be either — 100 Myr old clusters or clusters in the 
red loop of the models at ages around — 4 Myr. 

We notice that there are some objects with V — I > 1 
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Figure 6. Johnson-Cousins B — V versus V — I color-color di- 
agram (top) and V versus V — I color-magnitude diagram (bot- 
tom) for the objects in zone 5 (cluster candidates) that have been 
cross-correlated wi th the WFPC2 (B, V) photometry presented by 
lAloisi c t al. (1999). Full symbols denote objects located far away 
from the MB and SB which could be background galaxie s. Over- 
plott ed are the Padova simple stellar populations (IGirardi et al.l 
l2OT0t) for a Z=0.0004 metallicity, ages from 1 Myr to ~13 Gyr, 
and masses of 10 3 , 10 4 , and 10 5 M Q . Small dots indicate the 
Log(age)=6.7, 7, 8, 9 and 10 models. The red arrow represents 
an A v = 1 extinction vector (see text). 



which are redder than the oldest SSPs. We can ex- 
clude that they are background galaxies because they 
are all concentrated within IZwl8's MB, except for two 
of them, which are more "peripheral" . Another possi- 
bility is that these objects are highly reddened young 
embedded clusters; however, in this case they should 
have Ay > 3 (see reddening vector in Figure 0, and 
thus absolute magnitudes brighter than My ~ —10.75. 
Such high luminosities would imply a super star clus- 
ter (SSC) classification according to the definition by 
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Figure 7. Masses versus ages derived for the candidate star clus- 
ters with the Padova SSPs. The objects with V — I > 1 in Fig. [6] 
have not been included. 



iBillett et all ([2002). However, it is unlikely that IZwl8 
host so many SSCs, considering that NGC 1569, the 
dwarf with the highest known c urrent SFR, presen ts only 
4 of them (jHunter et all 120001 lOriglia et aLll200l . Sev- 
eral studies have discussed st ochastic effects on the inte- 
grated colors of clust ers (e.g. iSilva- Villa fc Larsenll2011l : 
iMaiz Apellanid 120091) . Stochasticity in the IMF sam- 
pling of a low-mass cluster can heavily change the ob- 
served colors with respect to predictions, and in particu- 
lar, red supergiants make appear the cluster colors red- 
der. iCervino fc Lur idiana (2004) investigated this issue 
and derived a luminosity limit under which stochastic ef- 
fects are important. Comparing the luminosities of our 
clusters redder than V — I ~ 1 with their predictions for 
ages older than 100 Myr, we found that our cluster lu- 
minosities are well above that limit. Therefore, the most 
likely explanation is that these objects are supergiants 
with ages younger than ~20 Myr, implying that the MB 
has been forming stars at very recent epochs. 

For the zone 5 objects with V — I < 1 we have inferred 
ages and masses using the Padova models. This was done 
minimizing the difference between model and data col- 
ors (B— V, V— I, and B— I) according to a x 2 criterion, 
varying the age between ~1 Myr and ~13 Gyr. We ac- 
counted for a Galactic reddening of E(B-V)=0.032, but 
did not solve for possible internal extinction, because of 
the small number of bands (3) and the relatively short 
wavelength ran ge covered by our data. Following the 
prescriptions of lAvnil (|1976| ). the 1 a errors were derived 
as the limits of the 0.68 confidence level region defined 
by X 2 < Xmin + 1- Cluster masses were obtained from 
the mass-to-light ratios of the best-fitting models in the 
I band. The results are presented in a mass versus age 
diagram in Fig. [7] Cluster ages span the whole Hubble 
time, even if we notice that the majority of clusters have 
ages younger than a few hundred Myr. Cluster masses 
range from 2 x 1O 3 M to 4 x 1O 5 M . The observed cor- 
relation between cluster ages and masses is an obvious 
effect of the magnitude limit adopted for zone 5. 
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We notice that our analysis provides 5 "old" (age>l 
Gyr) clusters, with masses between 6 x 10 4 and 4 x 
1O 5 M , and totaling ~ 10 6 M Q . This implies that ~3% 
of the stellar mass in IZw 18 is in the form of old clusters, 
which sounds unlikely given the fact that this ratio has 
been found to be considerably lower in other dwarf ir- 
regulars with evidence for significant SF at ages >1 Gyr 
(e.g., <0.3% in NGC 4449, Annibali et al. 2011). This 
suggests that the majority of the candidate old clusters 
in IZwl8 are instead resolved supergiants. 

Concluding, our analysis seems to suggest that the ma- 
jority of the objects in zone 5 are indeed individual stars, 
even if we can not exclude that some of them are unre- 
solved stars clusters. Nevertheless, even in this scenario, 
our analysis provides evidence for a particularly active 
SF over the past hundred Myr. 

5. DISCUSSION OF THE RESULTS 

The evolutionary properties derived above for the stars 
measured in IZwl8 are similar to those of other metal- 
poor starbursting dwarfs. For instance, once the differ- 
ence in distance modulus A(m — M)o ~ 0.7 is taken into 
account, magnitudes and colours of the various types of 
stars on the red plume of IZwl8 are in perfect agreement 
with those of the corresponding stars in SBS1415+437, 
the other very metal-poor BCD w hose stellar popula tions 
have been resolved by HST/ACS (jAloisi et al.ll2005h . In- 
deed, the CMD of SBS1415+437, whose intrinsic dis- 
tance modulus is 30.66, shows an unmistakable TRGB 
at / = 26.67 and the horizontal finger typical of carbon 
stars at / = 25.78 (jAloisi et all 120051) . SBS1415+437 
has a red plume more populated than IZwl8, as if it 
had been more active in previous epochs. In turn, the 
blue plume of IZwl8's MB is more prominent than that 
of SBS1415+437, suggesting that its current and recent 
star formation activity is proportionally higher. 

The results described in the previous sections allow us 
to draw a qualitative picture of the evolution of the MB 
and SB. Both components contain old and young stars in 
all evolutionary phases, but it is apparent that the two 
stellar systems are distinct from each other and have had 
different histories. The quantitative star formation his- 
tory (SFH) of the two bodies of IZwl8 will be derived 
fro m the CMDs using the synthetic CMD method (see 
e.g. iCignoni fc Tosil 120101 and reference s therein) and 
will b e presented in a forthcoming paper (|Annibali et al.l 

1ml). 

5.1. The Main Body 

The spatial distribution of stars in different evolution- 
ary phases displayed in Figures [2] and [5] suggests that 
the MB is indeed one single system and not "two galax- 
ies separated and interconnecte d by a narrow luminous 
bridge" as originally thought by Zwicky (196§). We can 
however understand how he was misled by the 8-shaped 
distribution of the young blue stars. 

Young massive stars (zones 1 and 3 of Figure [T]) cover 
the entire MB. Zone 3 stars are concentrated on the NW 
and SE star forming regions, following the traditional 
scheme of starburst dwarfs, where massive young stars 
usually group together and are concentrated in the inner 
regions, although not necessarily at the very center of 
the galaxy. The NW region contains a larger number of 



bright massive zone 3 objects than the SE region. This 
suggests a more recent/higher star formation activity in 
the NW clump. 

Zone 1 stars seem (see Figure [5]) intriguingly to avoid 
the areas where zone 3 stars are located. Is this be- 
cause they arc different objects? Since they are redder, 
zone 3 stars might be more evolved or more metal rich 
than zone 1 stars. We rather believe that either zone 1 
stars lie on gas filaments and are slightly blue-shifted by 
residual gas contamination increasing the F555W flux, 
or that differential reddening slightly affects zone 3 ob- 
jects. Throughout this paper we have only considered 
the foreground Galactic reddening, E(B — V) = 0.032, 
but IZwl8 does contain a lot of gas and presumably 
dust and molecule s, although in quan tities limited by 
its low metallicity (jCannon et al.ll2002l h It is thus plau- 
sible that absorbing interlopers are more abundant in 
the innermost, star forming regions and that they make 
the stars of those regions appear redder. In this frame- 
work, zone 3 stars are the same kind of objects as zone 
1 stars, only affected by internal extinction of the order 
of AE(V — I) ~0.1. 

Lower mass young stars (zones 2 and 4) could help 
in discriminating between the above interpretations, but 
they are faint and are much more severely affected by 
crowding than by any other effect. Indeed, although 
likely to exist also where massive young stars are found, 
they are only detectable where the latter don't outshine 
them. As discussed in the previous sections, in FlO's 
catalogue many of these objects lie on the gas filaments 
protruding from the NW spur. While the vast majority 
of these faint F10 detections on the filaments are likely 
to be gas peaks misinterpreted as stars, we do believe 
that some of the brighter ones, which appear as actual 
stellar sources when individually and visually inspected 
in the images, are ind eed very young stars o r unresolved 
star clusters (see also iChandar et al.l 120051) . This leads 
to speculate that star formation is occurring at least in 
the portions of the filaments closest to the MB. Unfortu- 
nately, we don't know whether these filaments are only 
made of ionized gas or also of neutral one, since the cur- 
rent resolution of HI data (F. Fraternali private commu- 
nication, Lelli et al. in preparation) doesn't allow one 
to resolve them. Neither do we know whether the fila- 
ments are infalling, outflowing or comoving with respect 
to the MB. From van Zee's et al. (1998, their figure 
6) comparison of the H a and HI distributions, the fila- 
ments appear to be a distinct feature of the ionized gas, 
but the HI resolution is definitely insufficient. It is thus 
difficult to assess if star formation is physically plausible 
there; certainly it is extremely unlikely that it occurs in 
hot, possibly outfl owing bubble s, such as many of IZwl8 
ionized filaments (|Martinl|19 96'). 

Also the massive objects of zone 5 are centrally con- 
centrated, as usually found in starburst dwarfs. Whether 
they are resolved supergiants or unresolved star clusters, 
they are the signature of strong star formation activ- 
ity. The evidence that they are more abundant in the 
NW than in the SE clump suggests that the activity 
in the former is higher than in the latter. The NW 
clump is also known to contain Wolf-Rayet stars ( see e.g. 
iHunter fc Thronsonl HIM llzotoy fc Thuaiil [1991 Yet 
from emission line spectroscopy. Ilzotov fc Thuan ([1998D 
suggested it to be not as young as the SE clump: 10 
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Myr the former and 5 Myr the latter (and 20 Myr the 
SB). The CMDs of Figure[4]do not seem to confirm their 
suggestion. Rather, our CMDs suggest that the star for- 
mation activity in the NW is more recent than in the 
SE. The circumstance that, of the variable stars reliably 
characterized by F10, the two massive ULP Cepheids VI 
and VI 5 are located in the NW region, while the lower 
mass Classical Cepheid V6 and the fainter LPV V7 are 
located in the SE region, favours the view that the NW 
on average contains younger stars than the SE. 

All the lower mass evolved stars are likely to be evenly 
distributed in the MB. RGB stars seem to be located 
more externally than AGBs, but this is most probably a 
selection effect. We ascribe their apparent confinement 
in the outer regions to the severe central crowding which 
prevents their detections. Similarly, we interpret their 
asymmetrical horseshoe distribution as due to the ex- 
tinction from the thick gas fi laments of the northern ar- 
eas fsee lCannon et all [2002. for the spatial distribution 
of dust). There is no doubt that faint objects are more 
easily detected in the SE external portions of the MB. 

The resulting overall scenario of the MB consists of a 
system with two active, distinct regions of current star 
formation, one (the NW one) stronger and younger than 
the other. These two active clumps are embedded in a 
fairly homogeneous body of stars of all ages, from a few 
Myr to several (possibly 13) Gyr. 

5.2. The Secondary Body 

The SB is in a rather different situation. As mentioned 
ab ove, the weaker emissio n lines of its spectrum already 
led llzotov fc Thuanl Jl998) to suggest that it is older than 
the MB em i ssion regions. WFPC2 photometry allowed 
lAloisi et al.l (|1999f) to quantitatively assess that its recent 
SFR is much lower than in the MB. The ACS photometry 
confirms these results. 

Our CMDs show that the SB does contain young blue 
plume stars, but they are more likely stars at least 10 
Myr old on the blue loop phase. Stars more massive 
than ~ 15M & do not seem to be still alive. The SB also 
appears to contain much less ionized gas than the MB: a 
further confirmation of its low recent activity. 

On the other hand, low and intermediate mass stars 
are at least as numerous in the SB as in the MB, with 
the objects of zone 4 (low mass MS stars) even outnum- 
bering those in the MB (see Table 1). This is a strong 
indication that at epochs earlier than ^100 Myr ago, the 
SB must have been actively forming stars as, or more 
than, the Ml£|. If we compare the CMD of the MB SE 
region (Figure |4} with that of the SB (Figured]), we see 
that the two red plumes are similar, both in the RGB 
and the AGB sequences. The blue plume, although con- 
taining similar numbers of stars, clearly shows that those 
in the SB are more evolved than those in the SE, that 
reach brighter magnitudes and bluer colours. The AGB 
and the RGB of the SB are well defined. The colour ex- 
tension of the latter, combined with the small metallicity 
range of the system, indicates that stars with the oldest 
possible ages (13 Gyr) are present. 

We conclude that the SB has started its star formation 

1 A deeper and more detailed analysis of the relative star forma- 
tion activity in both bodies, takin g; into acount completeness will 
be present in Annibali et al. (2011). 



activity up to ~ 13 Gyr ago, and has remained active 
since then, at rates comparable with those in the MB 
until relatively recently. However, in the most recent 
epochs (i.e. later than 10 Myr ago) its SFR has dropped 
to values significantly lower than in the MB. 

6. SUMMARY AND CONCLUSIONS 

We have studied the evolutionary properties and spa- 
tial distribution of the stars resolved in IZwl8 by 
HST/ACS proprietary and archival photometry. To this 
aim we have analysed both the A107 and the F10 cata- 
logues and found consistent results, independently of the 
different selection criteria adopted in the two approaches. 
Nonetheless, to be conservative, we have based our con- 
clusions on the A107 most restrictively selected catalogue. 

The comparison of the CMDs resulting from the ACS 
photometry with stellar evolution models (tracks and 
isochroncs) indicates that stars of all ages are present 
in the two bodies. In spite of the impossibility of reach- 
ing the oldest MS turn-off in a galaxy 19 Mpc away, 
the extremely low metallicity of IZwl8 relaxes the age- 
metallicity degeneracy and makes the red edge of the 
RGB a signature of the presence of stars up to possibly 
13 Gyr old. The oldest detectable stars are best visible 
in the SB and the SE portion of the MB, where crowding 
is less severe, but are present also in the rest of the MB, 
although measured with larger uncertainties. 

We have found that the stars are homogeneously dis- 
tributed over the two bodies, with the younger ones 
more centrally concentrated, as always found in starburst 
dwarfs, and the older /fainter ones spread out to the sys- 
tem periphery, also because of selection effects in their 
dctectability. From the maps of the stars spatial distri- 
bution (Figured]) we conclude that old and intermediate 
age stars are actually distributed homogeneously over the 
two bodies. 

Summarizing the results presented by A107 and F10 
and those obtained here, the main properties of IZwl8 
stellar populations are: 

• The distance independently estimated from the 
TRGB (A107) and the Cepheids (F10) is - 19 Mpc. 

• As far as stars are concerned, the SB is completely 
separated from the MB: no star is detected between 
the two bodies. If tidal interactions are at work, 
they don't appreciably affect the distribution of the 
stellar populations of the two components. 

• IZwl8 contains stars in all the evolutionary phases 
visible at its distance: massive and intermediate 
mass stars on the MS, blue and red supergiants, 
AGB and TP- AGB stars, carbon stars, RGB stars. 
Hence, it hosts stars of all ages, from a few Myr up 
to possibly 13 Gyr. 

• The MB hosts two bright star forming regions (NW 
and SE), embedded in a common and homogeneous 
environment of older stars. Some stars may have 
been recently formed on the gas filaments protrud- 
ing from the NW portion of the MB. 

• While old and intermediate age stars are present 
in similar amounts in the MB and in the SB, the 
youngest stars are only found in the MB and mostly 
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concentrated in its NW region. The SE region has 
also been active recently, but at a lower current 
rate; the SB has had definitely less star formation 
than the MB in the last tens of Myr. 

• Of the 34 candidate variable stars studied by F10, 
19 and 15 are in, or around, the MB and the SB 
respectively. All the confirmed ones are however in 
the MB: the two ULP Cepheids in the NW region; 
the Classical Cepheid and the two LPVs in the SE 
region. 

• At IZwl8's distance, star clusters are not distin- 
guishable from point-like sources. It is thus possi- 
ble that some of the brightest objects in the main 
body are actually clusters and not individual stars. 
Comparing their colors with simple stellar popu- 
lation models, we find that the candidate clusters 
would span ages from a few Myr to a Hubble time, 
with the majority of them being younger than ^200 
Myr. However, the fraction of stellar mass in old 
(age>l Gyr) clusters would be significantly higher 
than what found in other dwarfs, suggesting that 
the majority of the candidate old clusters are in- 
stead resolved stars. The reddest (V — I > 1) 
brightest objects cannot be unresolved clusters and 
must be supergiants with ages younger than ^20 
Myr. 

These results show that IZwl8 is not unique in the 
sky, but shares the evolutionary properties typical of 
BCDs and, more in general, of late- type dwarfs (see 
iTolstov et all 120091 for a review). It is not experiencing 
now its first burst of star formation, but has been forming 
stars over many Gyrs, most likely over the whole Hub- 
ble time. This implies that its extremely low metallicity 
cannot be explained by a low chemical enrichment due to 
extremely recent star formation. At best, one may hope 
tha t the SFR has been low until relatively recently (see 
e.g. lLegrand1l2000ft . although the conspicuous presence of 
stars of all ages does not make this likely. 

The same problem obviously concerns all the other 
metal-poor BCDs where s tars many Gyr old have been 
found, like SBS1415+437 (lAloisi et al.ll2005lh We recall 
that no genuinely young galaxy has been found yet by 
anybody. All the galaxies whose stellar content has been 
resolved so far contain stars as old as the lookback time 
allowed by the available instrumentation, wh atever their 
morphological type a nd metallicity (see e.g. iTosil 120091 : 
Cigno nTfc Tosil [201Ct and references therein). Hence, 
since stellar nucleosynthesis has had plenty of time to 
pollute them, the only ways to explain the low observed 
metallicity of these galaxies are either to dilute suffi- 
ciently their insterstellar medium by accretion of primor- 
dial or very metal poor gas, or to remove the stellar nucle- 
osynthesis products through galactic winds triggered by 
supernova explosions, or both. Chemical evolution mod- 
els of late- type dwarfs and BCDs have s hown long ago 
dMatteucci fc Tosil [19851 iPilvuginl [1991 iMarconi et all 
1 1994ft that infall of metal poor gas alone is not sufficient 
to reproduce the low metallicity of starburst dwarfs, and 
that galactic winds must be invoked too. 

Observational evidence for galactic w inds has been 
found for a handfu l of st arburst dwarfs (Boman s et al.l 
120071: iGrimes et all 12009ft . and IZwl8 is indeed one of 



those where the outflowing (Martin 1996) ionized gas as- 
sociated with the most recent SF activity may indeed 
have a speed higher than the escape velocity and be able 
to remove metals from the system. It would be crucial 
to get detailed information also on the cold gas, see how 
it is spatially distributed, and what its kinematics is. 

Understanding the chemical and dynamical evolution 
of IZwl8 represents a key step to approach the evolution 
of all dwarf galaxies as well as that of more massive sys- 
tems which may have formed from the assembly of small 
objects like it. 
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